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A method and apparatus for determination 
of heart rate, heart stroke volume, and cardiac out- 
put from thoracic bioimpedance signals and elec- 
trocardiograms. A unique bioimpedance electrode 
arrangement is employed, and the bioimpedance 
signals are corrected for gain phase-frequency dis- 
tortion through the use of sinusoidal test signals 
through the measuring or detection electrodes to 
identify distortions and correct for same during ac- 
tual measurements. Time-derivative bioimpedance 
signals are employed, the power spectrum cal- 
culated, and a novel autocon volution procedure 
used to emphasize the heart rate harmonic. Breath 
waves and other signals not indicative of the pa- 
tient's cardiocycles are removed. Left ventricular 
ejection time is derived from the bioimpedance sig- 
nals, and an improved version of Kubicek's equa- 
tion is employed to derive heart stroke volume and 
thus cardiac output. 
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NON-INVASIVE MONITORING OF HEMODYNAMIC PARAMETERS USING 

IMPEDANCE CARDIOGRAPHY 



5 TECHNICAL FIELD 

The present invention relates generally to cardiac monitoring and specifically 
to the determination of heart rate (HR), heart stroke volume (SV), and cardiac 
ou *put (CO) according to detection and complex analyses of thoracic bioimpedance 
and electrocardiograph (ECG) signals, which pennit precise detection of the start of 
10 left ventricular ejection. 

BACKGROUND ART 
Heart rate is the number of times the heart beats per minute. Heart stroke 
volume is the volume of blood pumped during each heart stroke. Cardiac output is 
the volume of blood pumped in one minute and is generally considered to be the 

15 most significant gauge of cardiac fitness. Physicians must frequently rely upon such 
cardiac parameters to diagnose heart disease, to assess a patient's overall health, to 
determine the most appropriate method of treatment, and to quickly discover sudden 
lapses in cardiac performance. 

The currently existing methods for measuring cardiac output and other 

20 cardiac parameters may be divided into two categories: invasive and noninvasive. 
The invasive methods require that a medical practitioner insert a measuring device 
into the patient's body, such as a catheter in the throat, and present numerous 
disadvantages to both patient and physician. The patient must often endure 
substantial pain and discomfort and the physician must perform a relatively 

25 complicated procedure and occasionally expose himself or herself to the risk of 

contact with infectious blood. The noninvasive methods currently in use represent a 
major advancement, but still have significant shortcomings. Most take 
measurements using ultrasound, phonocardiography, or electrical bioimpedance in 
order to calculate cardiac parameters. 

30 The methods which employ bioimpedance measurement involve placing a 

plurality of electrodes on a patient's skin (predominantly in the thoracic region), 
generating a high frequency, low amplitude electric current from certain of the 
electrodes into the patient's body, measuring the changes in the electrical impedance 



SUBSTITUTE SHEET (RULE 26) 



W ° 97/37591 PCT/SG97/00013 



-2- 



M=m, ^ over limCi ami ronehU]ig tte , n dKtttoi 

w»h cardiac parameters. I""*™* 
The manner of arranging te m teww^y ta 

~— * » various anatomical flctOT5 , ^ h 

«. of .he tad, ,„ ^ ^ ^ raeaMKd 

ml.rc^^carfhcpmn^. Many of the electrode configurations 
in use rail to adequately „, ^ fc ^ 

10 measurement. Moreover, a few electrode arrangements require the use of band 
dec**, e. g ., influ^ ^ dMrote ^ jj a]K) ^ deanjdes 

"* * """ v- BO. 1) These band elecuode, typically wrap ' 

«»nd a paier,, Uke , be„ and limit ^ , 0 ^ ffi 

condition during relation procedures. The mc-^ , BK i«ed 
Mb respiration aiso make „„, yKy taeon^ wta m ^ 

ttficfc and chest. 

Perhaps the most significant proNem Kit|l ^ 
bmrmpedance methods is the imprecise natdKMaOcal ttaivabon of oardiac 
1-rarnetnrs from bioimpedance measurements. The ventricniar ejection tf> (VET) 
h a measure of the toe between the opening and dosing of the aortic ^ 
*mng the sysmfc-diasn* cycle of the heartbeat and i, must be .alcnlated a, m 
■mediate step i„ defcrmUUng cardiac stxote volume. The prior „ doe, not 
«** a method for determimng vemrkutar Section time with sufficient 

Furthermore, the prior art fails to account for the fact that VET is not a single ' 
-*». 'Kere is actualiy a left Vffr and a right VET. It has been shown 

fat ft. rime-denvative impedance signai i, actually ptopo,,^ K tlle ^ ^ 

Mood flow ejected b, the left ventricle. The measurement of left VET and right 

VET for most patta* are genenllly nty ^ „„, ^ ^ 

-hem can create errors in bioimpedanc. readings under the methods currently in use 
Fnrtnermore, the classic algorithm for ejection Stan time is elaborate and 
weil oniy for healrhy patients at rest. „ i, M accl , rate fot ^ ^ 

Pbys.c.1 training or other physic, sues, or for critoUy i„ pafe,,,. mi M lhosc 

typically in intensive care units. 



WO 97/37591 PCT/SG97/00013 

-3- 

The conventional equation for deriving stroke volume from bioimpedance 
signals has become known as the Kubicek equation and is given as follows: 



SV = ROIZJ 1 ■ AZ 

5 where SV is heart stroke volume, R is blood resistivity, L is the distance between 
the inner and outer voltage sensing electrodes, Zo is the mean thoracic impedance 
determined from the inner voltage sensing electrodes, and AZ is the impedance 
change due to blood influx. Kubicek's estimation of this value is 

AZ = (VEDidZfdt)^ 

10 where VET is the combined left and right ventricular ejection time, and (dZ/dt^ is 
the maximum negative slope change of the time-differentiated impedance signal. 
Most bioimpedance cardiac monitoring systems use some form of the Kubicek 
equation. 

Without further refinement, however, the Kubicek equation frequently gives 
15 inaccurate measurements. This is due in part to the fact that both ventricles 

contribute to impedance changes, and so Kubicek's calculated ejection time (VET) 
cannot be associated with a particular, specifically the left, critical ventricle. 
Concurrently, Kubicek's AZ estimation becomes invalid when strong left-right 
ventricles asynchronism is observed. As a result, Kubicek's SV calculation is often 
20 proportional to, but not equal to, the actual heart stroke volume and must therefore 
be multiplied by some correlating constant. In addition, the prior art does not 
disclose a method for adjusting R in accordance with the fluctuation of a patient's 
hematocrit (red blood cell count). The adjustment of R is especially important in 
patients undergoing blood infusion. 
15 Many of the methods for bioimpedance cardiography require that the patient 

hold his or her breath during each measurement because respiration causes 
interference in (he bioimpedance signal. Such methods are inconvenient for some 
patient? and completely useless for odier patients who are unconscious or otherwise 
unable to hold their breath. Some of the more recent methods include signal 
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DISCLOSURE OF INVENTION 
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used to measure the bioimpedance and ECG. The computer system communicates a 
sinusoidal test cum to the transducer receiver and may then measure and record the 
gain-phase-frequency distortions created by the receiver. The "real' bioimpedance 
signals subsequently received by the computer system through that transducer may 
5 pass through filters which remove the transducer's characteristic distortions as weD 
as bieath and movement artifacts. Accuracy of gain-frequency response and phase- 
frequency response may be corrected to within 5%. 

The invention may further comprise a method of signal processing the ECG 
signal to determine QRS complexes (characteristic heartbeat waveforms) and check 
10 point positions for use in refining the bioimpedance signal calculations. The signal 
processing may involve sampling the measured ECG signal and filtering it to 
highlight the positions of the QRS complexes. Peak-to-peak amplitudes may be 
recorded and a threshold amplitude may be calculated. Defect-free QRS complexes 
are selected using the threshold amplitude. Finally, additional analysis of selected 
1 5 events may be performed to determine check point positions. 

The computer system of the invention may derive and save in memory a 
time-derivative bioimpedance signal. The curve generated by the lime-derivative 
bioimpedance signal may be platted with respect to time and represents repeated 
cardiocycles. The power spectrum of the bioimpedance signal may be calculated 
with discrete Fourier transforms and studied to estimate the patient's heart rate and 
to identify the fronts of each cardiocycle. A novel mathematical autoconvolution 
procedure may be used to emphasize the heart rate harmonic in the time-derivative 
bioimpedance signal. 

Breath waves may be removed by generating an envelope within the power 
spectrum in which cardiocycle signals should be found and then removing those 
signals which lie outside the envelope. 

The computer system may employ a new method of deriving the effective 
left ventricular ejection time (ELVET) from bioimpedance signals and ECG signals. 
In particular, precise detection of the ejection start time is based on complex 
analyses of the bioimpedance and ECG signals. The derivation of ELVET is a 
detailed, multi-step analysis which involves finding various points on the time- 
derivative bioimtKdance curve based on corresponding points on the ECG curve and 
determining which of these point, most accurately reflect cardiac events. It requires 
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BRIEF DESCRIPTION OF DRAWINGS 

20 patient;^' 1 " ' * * * C,n * on a 

FIGS. 2A, 2B, and2C statically depict spot-type electrode placements 
according to the present invention; 

^ 3A ana ^al zones of interest for placement of measuring 
chest electrodes according to the present invention; 

^ patient ^ ^ *" P ° Siti0nine ° f Mat ° micaI 20065 ° f HG - 3 on a 

HG. 4 illustrates the structure of an individual chest electrode according to 
the present invention; 

FIG. 5A is a gnph of (ho jain »„„ frcq^^ characteristic, of an A-filte, 
accort.06 to .he Beaton c^ 1( ^ ,„, ^ ojllKM|i , |ioM >niJ 

suppression; 

HG. SB is a jrap,, of „« ^ vtnw of a 
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FIG. 5C is a graph of filtration employed to remove breath waves' power 
spectrum from the summary (breath and cardiowave) power spectrum in the present 
invention; 

FIG. 6A is a branch chart of the methodology of the present invention; 
5 HG. 6B is a flow chart of signal processing and depicting the apparatus 

employed in performing the method of the present invention; 

FIG. 7A depicts the gain-frequency characteristic G(f) of a bioimpedance 
transducer according to the present invention; 

FIG. 7D depicts the phase frequency characteristic P(f) of a bioimpedance 
10 transducer of the present invention; 

FIG. 7C depicts the gain frequency characteristic * of an "R" or 

Gif) 

restoring filter corresponding to the characteristic of FIG. 7A as used for GPF 
corrections according to the instant invention; 

FIG. 7D depicts the phase frequency characteristic -P(f) of an M R H or 
15 restoring filter corresponding to the characteristic of FIG. 7B as used for GFF 
corrections according to the instant invention; 

FIG. 8A depicts the resulting gain frequency characteristic G(f) of a 
transducer signal as corrected by an R-filter according to the present invention; 

FIG. 8B depicts the resulting phase frequency characteristic P(f) of a 
20 transducer signal as corrected by an R-filter according to the present invention; 

FIG. 9 is an illustration of a time-derivative bioimpedance signal generated 
according to the present invention, plotted with respect to time and depicting 
exemplary impedance changes due to blood influx (AZ) over the same time period; 

FIG. 10 is a graph of the electric potential traversing the heart of a patient 
25 plotted with respect to time (electrocardiogram), depicting die peak-to-peak 
amplitudes (E, - front, E, - back) of the QRS complex in a single heartbeat; 

FIG. 11 is a scatter diagram of all peaks measured in an electrocardiogram 
over a period of ten seconds. The coordinates of the graph are (E,, Ej, and the 
QRS complexes are within the dotted circle in the lower right hand comer, 
30 FIG. 12 is a graph of the ECG signal and the time-differentiated 

bioimpedance signal plotted against time indicating check-points P, S„ and S b on the 
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ECO signal, relative to check-points Q, S, and A on the time-differentiated 
bioimpedance signal; 

HG. 13 is a graph of ECG) bioimpedance m ^ time ^ ffcrcntiated 
b^pedance <W dt , ^als with respect to time Moating the bioimpedance 
difference between the S and Q points. Z„ in a healthy patient; 

HG. 14 ls a graph of ECG, bioimpedance Zft), and time-differentiated 
bm,mpeda„ce signals with respect to time indicating the bioimpedance 
difference between the S and Q points, Z„ in a ^ with ^ ^ ^ 

WG. 15 is a graphic depiction of the correlation of the inventive method of 
detomining ELVET to an ideal, ultrasound deterrnination in comparison to the 
classic Knbicek algorithm methodology. 

BEST MODES FOR CARRYING OUT THE INVENTION 
The first step in the present invention involves taking bioimpedance 
measurements over segments of tissue on a patient's body. Electrodes must be 
Placed at appropriate points on the surface of the skin to generate a high frequency 
low amplitude electric current and to detect changes in the generated current after ,'t 
passes through the segments of tissue (see FIGS. 2A, 2B. 3 and 4). The electrodes 
are "spot electrodes" rather than "band electrodes ' in order to maximize the free 
area on the patient's body. The spot electees are preferably of the disposable, 
one-use type. The patient thus has increased freedom of movement and medical 
practmoners have more access to the patient's skin for other medical procures 
such as the introduction of catheters and the administration of anesthesia 
25 The bioimpedance electrode system employs a total of six electrodes: a pair 

of detecting Measuring) electrodes 20 at the xiphoid process level, a pair of 
detectmg (measuring) electrodes 22 positioned laterally on the neck, an influencing 
electrode 24 on the left leg, and an influencing electrode 26 on the forehead 

The mfluencing electrodes 24, 26 may be standard ECG spot electrodes with 
30 a contact area of 2 centimeters by 2 centimeters (cm). The upper influencing 
electrode 26 is preferably placed on the middle of the forehead, at the mid-line 
thereof. The lower influencing electrode 24 is preferably placed on the left knee or 
somewhere below the left knee such as the left foot. If necessary, the lower 
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influencing electrode 24 may also be placed above the knee level, provided that the 
following condition is satisfied: 



L > SR. 



where L is the distance between influencing electrodes and R is the radius of the 
chest. Tiie left leg is used instead of the right leg to account for the anatomic 
asymmetry of die heart. The physiological positioning of the aortic arch, through 
which a significant concentration of electrical potential power lines pass, makes the 
left leg most suitable for the lower influencing electrode. The arrangement of 
influencing electrodes in this manner guarantees the uniform distribution of 
influencing current power lines between measuring electrodes and thus helps to 
minimize the error in the final cardiac parameter measurements. 

The upper pair of measuring electrodes may also be standard ECG spot 
electrodes with a contact area of 2 cm by 2 cm. These electrodes are placed 
symmetrically along the lateral lines of the patient's neck about the perimeter of the 
patient's neck 27, a distance S above the base of the neck. The distance S is 
defined as the distance between the base of the neck and the center 302 of electrode 
22 and is preferably approximately 4 cm. The base of the neck is defined to be 
located at the point of maximum curvature of the lateral lines of the neck. Placing 
the upper measuring electrodes 22 in this area avoids the error that would otherwise 
result from the noolinearity of electrical power lines at the neck-chest junction. 

The lower chest pair of measuring electrodes 20 each have a contact area of 
12 cm 2 to 30 cm 2 . If this contact area is either reduced or enlarged, the heart 
stroke volume will be underestimated A contact area of less than 12 cm 2 provides 
insufficient depth of measurement, a particularly serious problem with larger 
patients, and a contact area of more than 30 cm 1 causes the measurement to extend 
into additional anatomical regions. 

The individual chest electrodes 20 are each preferably comprised of a set of 
four standard ECG spot electrodes 28, each with a contact area 29 of 2 cm by 2 cm 
with the top pair 34 of spot electrodes 28 of each chest electrode at xiphoid process 
level 38 (see FIGS. 2B and 4). All contact areas 29 are connected wiLh foil or wire 
31. The distance G separating adjacent spot electrodes is approximately 5 cm. The 
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featutes of zones I and m of die bod, (see PIOS. 3A and 38, and ins„ ra , 
munmum of error rince the .n^^, jj taken at an ademnKe depth a*l allow, 
5 ^^^^-^^oftSHe^pa^. TVCestc^oZ 

2B, and 2C) at the xiphoid process level 38. As shown in PIGS. 3A and 3B ' 
is » rones f a^ rv, wUhin an a« a^eriorly „ „ 

Tie toflnMctag electrode, 24, 26 genenr* , Ugh frequency, !o„ anrptode 
current into tho patient's body and the detecting electrodes 20, 22 measure the 

""^"^"♦^ ^^ctricaJtaped^ofthcB^ 
'"^^"■^dlfferer^^ 

measured cur™, . Tl« ek^, taf ^ of «« ^ ^ u a 
result of Wood flow, resr*Mion, and other factors. 

The present irrvemion also uses BCO signals concurrently measured wirj, the 
**l W ,he exception of post-mr^rern,*, signlil 
remove hardware artifacts described below, the ECG signal measurement is 
performed in a con^iona. manner, and so will „» be further described 

Apretal^uys^inthe present invention is to determta.tbegaiu-pn.se. 
««>«ency (GPF) chancterisfe of th. analog tornr. ^ mo ^ 
rapeuance converter) for subseonency use In signal processing. 

A computer syaem, which may comprise a specifically programmed general 
purpose compute, ^ a! a ^ ^ tfB ^ 

measurements from the detecting cleorode, 20, 22, define, bo,h th. impedance 
of te tnterceding ««« (bioimpedance) and the electric potential traversing th 
hew (ECG) as a faction 0 f time, and ultimaKly cremates HR (hear, ran,) Sv 
*7« — volume), and CO (Cardiac Output). The equation for draining SV 
30 « given as follows: 
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where K is a novel scale factor, P is the specific blood resistivity, L is the distance 
between the voltage measuring or sensing electrodes 20 and 22, 7^ is the mean or 
base thoracic impedance (determined from sensing electrodes 20 and 22) and AZ, 
the impedance changes due to blood influx. AZ is calculated as 

AZ - ELVET^dZfdO^Z^ 

where ELVET is the effective left ventricular ejection time, (dZ/dt)^ is the 
maximum absolute value of the thne-differentiated impedance signal obtained from 
the two measuring electrodes, and Z,., is a novel correction factor that lakes into 
account left-right ventricles a synchronism, equal to the bioimpedance difference 
between S and Q points (see FIGS. 13 and 14 and accompanying text). 

The K term accounts for variations in body compositions of different 
patients. In order to obtain a value for the K factor, the medical practitioner first 
measures the height and weight of the patient and the perimeter of the patient's neck 
27 and chest 36 using conventional means or commercially available ultrasonic 
measuring means. The medical practitioner then inputs these values into the 
computet system, which in turn uses the values to compute the effective cross- 
sectional area of the chest and the K factor. The effective cross-sectional area of 
the chest (SCHEST) is given by: 

SCHEST = (PCHESt 1 * (PNECK • POTEST) + PNECK 1 ) f 12n, 

where PCHEST is the perimeter of the patient's chest 36 and PNECK is the 
perimeter of the patient's neck 27. Then the K factor is calculated as follows: 



* " *o " *i ' (/SCHEST/ (H K » • W K >)), 
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where H is the patient's height, W is the patient's weight, and K^, K, , K„ K, are 
gender and age dependent and lie in ranges: 

*o « D-4]; JT, e [3-16]; £, e p>-l]; IT, e [0.1-2] 

The present invention can thus be used on patients of varying body constitutions 
5 without a loss of accuracy. 

The electrical resistivity of human blood is not a constant It varies among 
different individuals and even in the same individual at different times. Blood 
resistivity is particularly susceptible to fluctuation in patients undergoing blood 
infusion. As a result, an accurate system of bioimpedance cardiography must 
include means for continuously modifying the blood resistivity term of the Kubicek 
equation. 

A patient's specific blood resistivity depends largely upon his or her 
hematocrit. The relationship between these two values for capillary blood is as 
follows : 

P-I3J5+429-B, 

15 

where P is the specific blood resistivity and H is the capillary hematocrit This 
relationship is adopted from V.I. Arinchin et al., "Taking into account electrical 
resistance of blood will increase accuracy of chest tetrapolar rheography method," 
Journal of Pediatrics (U.S.S.R.) 1987, v. 7, pp.59-52. The hematocrit can be ' 

20 measured using any commercially available method. It may either be inputted into 
the computer system by the medical technician, or sent by electronic means directly 
from the hematocrit measuring device. 

The invention utilizes a novel method for processing the ECG signal, (after 
hardware artifact removal) comprised of the following steps: 

25 ® sign*! approximation from sampling points, 

GO special filtering to highlight the positions of QRS complexes, 

(iii) measuring peak-to-peak amplitudes for a given, recorded time 
interval, 

(iv) calculation of amplitude threshold, 
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(v) QRS selection with the calculated amplitude threshold, 

(vi) additional analysis of selected events lo determine check point 
positions. 

Step (i) above is desirable to increase the accuracy and reliability of QRS 
determination unless sampling frequency is extremely high. As a matter of 
practicability, such high frequency sampling is undesirable as consuming excessive 
processing time and memory capacity, and is impractical to effectuate the signals of 
restricted power spectrum. It is contemplated that two approaches to approximation 
may be suitable for use with the invention. First, it is well known that any signal 
s(t) with finite spectrum (and defining the highest harmonic as o> m =2rfJ, is fully 
described by its samples at points s(n*AT), where AT£l/24 is the sampling 
period, and n is an integer A precise approximation of such a signal is given by 
the equation; 



The same result can be obtained in another way by first calculating a Fourier 
transform of the signal s(t) and adding a small phase shift A?, to all harmonics, so 
that 



Then, after calculation of the inverse Fourier transform, the approximated values of 
s(t) at the points shifted by t from the original samples arc obtained. This latter 
approach is more efficient in calculation. Both methods permit lower sampling 
rates that lead to less consumption of memory and can be employed on demand to 
calculate a precise approximation of the original signal. 





Acp/2rc/ = t - const, Vf. 
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The next step is filtering the ECG signal to highlight the positions of QRS 
complexes, FIG. 10 illustrates a typical QRS complex, a signal peak with greatest 
amplitude measured from peak-to-peak in a single heartbeat. A symmetrical finite 
impulse respond (FIR) digital filter is calculated from the desired gain-frequency 
5 characteristic (FGC^). The desired GFC Btel is elaborated from analysis of power 
spectrum of QRS complexes and has a passband from 6 Hz to 22 Hz with the 
maximum at 12.5 Hz. Using a discrete Fourier transform, the desired GFC^ is 
converted to finite impulse characteristic according to an algorithm for filter 
synthesis described in V.S. Gutnikov, "Filtration of measured signals,* Leningrad 
10 Energoatomizdat (USSR) 1990, pp. 172-181, incorporated herein by reference. 

This filter passes the QRS complexes and suppresses breath and movement artifacts 
in the ECG signal, also as P and T waves. 

The next step is to calculate the peak-to-peak amplitude threshold and select 
valid QRS complexes. The computer system measures each local peak of the 
15 filtered ECG signal by its front (E,) and back (E,) amplitude fronts, see FIG. 10. 
For each local peak, E, is the measured from the peak's anterior, or leading local 
minimum to the. next nearest maximum, and E, is measured from the peak's 
maximum to its posterior, or trailing, minimum. FIG. 1 1 depicts a distribution 
(scatter diagram) of peaks by their (E,, E,) coordinates for a time interval or period 
of 10 seconds. FIG. U also shows the QRS complexes highlighted within a dotted 
circle in the lower right hand corner of the figure. Each peak is characterized by 
its (E,, E0 vector and amplitude A ip where 



20 



25 



\ = «V + <V> 

The computer system then searches the sorted {AJ array for the maximum 
difference between A, and A ltl . if the maximum is found for an exemplary *-th 
element, then (he amplitude threshold is calculated as T = (At+A^/2. Thus, the 
QRS complex is detected in point j if Aj exceeds the threshold T. The threshold T 
is adapted for each 10 second block or interval of ECG data as T. = T*., + <*T», 
where T, M is the adapted threshold for the previous data block, T» is the calculated 
threshold for the current block, and a is a parameter of adaptation in the range of 
30 0-1. 
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Each QRS complex identified by the use of the above "threshold" 
methodology is then farther analyzed within an interval of -50 to +200 milliseconds 
(ms) from the determined QRS position. For each QRS complex identified, the 
computer system determines the amplitude, sequence of peaks, and derivative of 
peak fronts to arrange 3 check points (FIG. 12): the start of the QRS complex (P- 
point), the maximum deviation from the base-line (S.-point. which coincides with 
peak R in a normal ECG), and the end of the QRS complex (Vpoint). These 
check points arc used to refine the analysis of the bioimpedance signal described 
below. 

The invention utilizes a novel method for processing the bioimpedance signal 
comprised of the following steps; 

(i) digital filtration and phase correction, 
(u) heart rate estimation, 
(iii) suppression of breath waves, 
15 (iv) determination of cardiocycles, 

(v) arrangement of check points, and 

(vi) selection of cycles without interference artifacts. 

The first part of the electronic filtration involves passing the signal through a 
"restoring" R-filter to achieve gain-phase-frequency (OPF) correction. The R-filter 
20 compensates for the distortions caused by the particular electronic transducer that is 
used to measure bioimpedance changes. It is well known that GPF characteristics 
of a bioimpedance transducer (FIG. 7) may greatly influence the shape of the 
bioimpedance curve. These interferences must be removed from the signal. The R- 
filter uses posterior signal processing to correct linear GPF distortions. It is 
constructed in such a way that the system of the bioimpedance transducer plus R- 
filter has GPF characteristics with zero phase shift and constant gain at the given 
range of frequencies, for example from 0.3 Hz to 30 Hz for the bioimpedance 
signal (see FIG. 8). Thus, correlation of the outputs of different bioimpedance 
devices may be achieved. With the exception of the filter parameters and operating 
30 characteristics specified as desirable or critical for the R-filter and other filters 
described herein, construction of same is conventional and within the ability of 
those skilled in the art, and so will not be further described. 



25 
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The first step of the GPF correction involves connecting the bioimpedance 
transducer with a source of an electricaUy-generated sinusoidal impedance signal 
and then measuring the output from the transducer. The electronically-generated 
sinusoidal impedance signal has an amplitude of 0. 1 Ohm to 0.2 Ohm with respect 
5 to a baseline, for example, 100 Ohm to 200 Ohm. Such a signal has been 

developed using a voltage-to-impedance converter consisting of a photoresistor, a 
photoemitter (photo diode), a power supply, and an analog-digital-analog (ADA) 
computer interface. The paired photoresistor and photoemitter are coupled inside a 
bght-protected housing so that the photoresistor changes its impedance according to 
10 the photoemitter's light intensity. The ADA conversion process includes digital-to- 
analog conversion of the mathematically-modded sinusoid with a frequency of 19 
kHz and analog-to-digital conversion with a frequency of 100 Hz, with 12 bit 
resolution. Through the interface, the computer produces a set of test sinusoidal 
signals with frequencies in the range from 0 Hz to 75 Hz and records the responses 
15 of the transducer. The operat.ng characteristics of the voltage-to-impedance 

converter indut!e ^ ^ sigMl 0 f 0 V to 5 V, an output signal of -0.1 Ohm to 0 I 
Ohm with a baseline of 100 Ohm to 200 Ohm (as previously noted). The OPF 
characteristic H(f) of the transducer may then be calculated from the spectrins of 
the initial test signals and the resulting responses of the transducer and presented as 
20 a graph or stored in an ASCII or other memory file. The system uses the calculated 
GPF characteristic H(l) of the transducer to calculate the "restoring- R-f,lter. The 
GPF characteristic of such R-filter can be formally written as 1/H(f) in a certain 
frequency range. The R-fH ter a i so provi ^ ^ ^ ^ 

frequency filters to provide suppression of random low and high frequency 
25 interference (see FIG. 7). Hie filtration with the R-filter may be done in the 
frequency domain using a Fourier transform. It is preferable to use a Gaussian 
window with Fourier transform to eliminate the boundary effects of the recorded 
signal: 

G(t) * eap[-2(ar/(2T)) *] 



30 



where 2T is a duration of recorded signal, t < T is a time, and "a" is a predefined 
constant preferably i n the range of 2.5-3. Multiplication of the Fourier image of the 
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rccordecJ signal by the R-filter's GPF characteristics results in suppression of GPF 
distortions and additional filtration of the signal. Reverse Fourier transformation 
and division by the Gaussian window may also be employed. The same steps may 
be employed in the time domain without Fourier transformation. The signal after 
5 R-filtxation is referred to as the "restored" signal, and this signal is used for further 
calculations. 

It should also be noted at this time that the identified GPF characteristics of 
the ECG connection are processed to remove hardware artifacts from the ECG 
signal in a similar manner to that described above for the bioimpedance signal. 
10 GPF correction of both the ECG and bioimpedance signals promotes true 
correspondence of time intervals and event times between the two signals. 

The next step in the bioimpedance signal processing is heart rate (HR) 
estimation. The present invention uses two ways to calculate HR. The ordinal way 
is to detect R-peaks on the ECG signal as described above and calculate R-H 
15 interval. The inverse value multiplied by 60 corresponds to heart rate. If ECG 
signal cannot be processed to detect R-peaks for some reason, the second way is 
used. In the second way, the power spectrum of the "restored* bioimpedance signal 
is calculated with discrete Fourier transform and used to estimate the patient's heart 
rate (HR). Very often a breath harmonic is the biggest one in the power spectrum 
20 of the bioimpedance signal. Consequently, it must be suppressed and the HR 

frequency response highlighted. A special transformation is used for this purpose. 
First, the power spectrum (PS) of the "restored" signal is multiplied with the gain- 
frequency characteristic of the A-filter (see FIG. 5A). This fdter differentiates the 
signal and additionally suppresses harmonics below a certain frequency preferably 
25 selected at a range from 1 Hz to 3 Hz because breath wave harmonics commonly lie 
below 2 Hz and HR harmonics above 0.8 Hz. The power spectrum of clear cardio 
signal consists of repeated peaks at frequencies HR, 2*HR, 3*HR, etc. 
Consequently, the following autoconvolution of the power spectrum will emphasize 
the heart rate harmonic: 

AS! (ft * PSa(J) PSa{2fl*PSa&fr.. 
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where AS is a result of the autoconvolution of the power spectrum and PSa(i) is 
the power of a K iven spectral line with frequency i which previously passed though 
the A-filter. It is preferable to have in the above product only the first three 
element? PSa, because the higher the frequency is, the lower the signal-to-noise 
5 ratio. The computer system uses the autoconvolution to search for Mas, the 
maximum value of ASl(i) in the range of 0.6 Hz to 5 Hz. The frequency 
associated with Mas is regarded as an estimation of HR. The estimation of HR is 
then used in additional filtration (see FIG. 5C) and caidiocyclc-identification 
procedures. 

10 The aberrations in the bioimpedance signal caused by respiration must be 

removed to increase accuracy and to insure proper identification of the cardiocydes. 
Usually, breath frequency is less than heart rate frequency, but breath waves create 
a power spectrum which overlaps the lowest harmonics of the power spectrum 
created by the cardiocydes. So, it is impossible to remove the breath waves' power 

15 spectrum from the summary power spectrum entirely (see FIG. 5C). Cardiac 
strokes are a more stable,* repetitive process in comparison with breath. 
Consequently, we can consider that their power spectrum consists of several narrow 
peaks. All power spectrum harmonics between main cardio waves' spectrum lines 
are combination of lateral slopes of these main spectrum lines and a noise power 

20 spectrum. If these internal harmonics are decreased, noise is mainly suppressed, 
and cardiowaves only slightly. As subsequently noted, these ideas form the basis 
for the algorithm for breath wave filtration. After estimation of the heart rate (HR). 
as described above, the first and second harmonics of the cardiowaves' spectrum arc 
determined. The local minimums at the power spectrum nearest to these spectrum 

25 peaks can be considered as their bounds. All harmonics below the lower bound of 
the second peak except for those within the first peak bounds are multiplied by a 
predetermined value less than 1 (as schematically shown in FIG. 5C). This results 
in elimination of breath wave amplitude, because the latter's power spectrum lias in 
the multiplied zone, but only slightly affects the cardiowaves. 

30 The next step in bioimpedance signal processing is cardiocycle detection. 

The invention also uses two ways to do this. The first way is to place the 
cardiocycle borders according to QRS complexes positions in die ECG signal. If 
the ECG cannot be processed, the second way is used. Jn the second way, the 
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" restored" signal as mentioned above is passed through the B-filter (FIG. 5B) to 
bighliglit the fronts of the cardiocycles. This is a differentiating filter with sinusoid- 
like frequency bounds. The pass band of the B-filter is adapted for frequency 
harmonics that produce a main contribution to cardiocycle front (preferably from 1 
Hz to 6 Hz). The area between a local minimum and the next local maximum at 
the signal passed through the B-filter Is regarded as a cycle front and described with 
a peak-to-peak change in time and a peak-to-peak change in amplitude. The 
computer system then generates a time-amplitude envelope by analyzing the first 
five (5) harmonics of the power spectrum generated by the signal after it passes 
through the B-filter. The cycle fronts are examined and certain ones are marked for 
further analysis if their peak-to-peak changes in time and peak-to-peak changes in 
amplitude are within the time-amphtude envelope. To increase the reliability of 
cardiocycles recognition, the computer system calculates mean and variance of 
peak-to-peak amplitudes for the selected regions. 

If the variance-per-mean ratio is less than a predetermined value, preferably 
0.3, then all of the marked fronts are transmitted to the next stage which involves 
the arrangement of check points and the selection of defect free cardiocycles. In the 
alternative, if the variance-per-mean ratio is greater than a predetermined value, 
then additional analysis must be performed. The additional analysis comprises the 
20 following steps: (i) the regions under examination are separated into two groups 
according to their peak-to-peak amplitude: those above the mean value ("the upper 
group") and those below the mean value ("the lower group"); (ii) the means, Ml 
and M2, and variances or standard deviations, VI and V2, are calculated for each 
group; (iii) the appropriate values are inserted into the following inequality for each 
25 group: 

M2+a>V2<Ml-<i'Vl> 



15 



where a is a predetermined value, preferably 1.96, and VI and V2 are standard 
deviations for the lower and upper groups respectively; and (iv) if the inequality 
holds true, then the regions in the upper group are taken to be the fronts of the 
30 cardiocycles ami the regions in the lower group are eliminated from further 
consideration, otherwise all selected regions proceed to the next stage. 
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Thc computer system identifies certain check poinfs in the time-differentiated 
bioimpedance signal to calculate the effective left ventricular ejection time, ELVET, 
as a preliminary step in determining heart stroke volume (SV). 

The present invention uses ELVET. a term which represents only the left 
5 ventricular ejection time, rather than Kubicek's VET, a term which represents the 
combined left and right ventricular ejection times. It is known that the value of the 
time-differentiated impedance signal is proportional to peak aortic blood flow 
ejected by the left ventricle. Therefore, the most accurate calculation of SV 
requires that LVET be used. LVET is calculated from the following equation: 

LVET=ELVET+LVPT, 

where ELVET is the time between the moment the left ventricular valve opens 
(S-point) and the moment that it begins to close (T-point), and LVPT is the 
protodiastoly time (the time it takes for the left ventricular valve to close). LVPT is 
not readily detectable with hemodynamic monitoring means because the changes in 
blood flow are insignificant during the protodiastoly time. For this reason, the 
present invention uses the product of ELVET and (dZ/dt)^, the maximum absolute 
value of the first chest impedance derivative with respect to time, to account for the 
lack of an LVPT measurement. 

The calculation of ELVET requires an analysis of the curve generated by the 
graph of Y(x), time-derivative bioimpedance, plotted with respect to time (see FIG. 
9). The computer system first finds a global maximum of time-derivative 
rapedance, Y(x), over a given cardiocycle and designates it as point A. The 
computer system then traces back in time from the A-point to the point in time 
corresponding to the S.-point on the ECG signal, a „d looks for abnormalities in the 
bioimpedance signal between those two points (see FIG. 12). The abnormalities of 
interest are: (1) dZ/dt zero crossing, (2) local minimum in dZ/dt, and (3) local 
maximum in the third derivative of the bioimpedance signal, dW. If there is no 
abnormality in the bioimpedance signal found in the time interval between points A 
and S., the cycle is considered defective and rejected from further consideration. If 
any abnormalities are found to the right of S>, the abnormality closest to S„ 
approached from the right is selected as the ejection start lime, S. Otherwise, the 
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abnormality closest to S h approached from the left is selected as the ejection start 
time, S. The usage of the ECO signal increases stability of the S recognition in 
complex cases. 

To identify the end of ELVET, point T, the computer system first finds 
5 point T 0 , which is taken to be either the first or the second local minimum after 
point A at the time differentiated bioimpedance signal, labeled as T, and T 3 
respectively (sec FIG. 9). The computer system chooses between T x and T 2 after an 
analysis of the depth (amplitude) of the curve at each point. If depth of the second 
minimum is greater than a predetermined fractional value of the depth of the first 
10 minimum, then T, is selected as T 0 . Otherwise, T x is used. To increase the 

stability of T„ detection in noisy signals, the invention looks for the back or trailing 
edge of the T-wave in the ECG signal. The back edge of the T-wave is detected 
from the local maximum at ECG signal next to QRS complex up to the next local 
maximum at the graph of curvature of ECG signal versus time. If one of the points 

15 Tj or T 2 is out of bounds of the T-wave's back, the other point is used as T 0 

regardless of its amplitude. The "effective end" of the ventricular ejection, point T, 
is then identified as the nearest local minimum before point T 0 on the graph of the 
curve generated by the second derivative of Y(x). ELVET is calculated as the time 
distance between points S and T. 

20 Hie inventive method of determining ELVET has been correlated to an 

ultrasound ELVET determination with a correlation coefficient of r = 0.86. 
Kubicek's classic algorithm gives only a correlation of r = 0.71. See FIG. IS for a 
graphic depiction of the correlations of the classic versus the new methodology to 
the ideal. Consequently, and in contrast to Kubicek, left ventricular ejection time as 

25 measured by the inventive method is measured in substantial isolation from the right 
ventricular ejection time. 

Normally, the point of the ejection start r S, coincides with the zero crossing 
of the tune-differentiated (dZ/dt) bioimpedance signal (sec FIG. 13), but this is not 
the case for seriously ill patients (see FIG. 14). For such patients the point of the 

30 ejection start is often placed at the "stair" or "abnormality" of the first front of the 
bioimpedance signal. This "stair" or "abnormality* characteristic (also leaned a 
•prewave") of an ill patient is commonly referred to as left-right ventricular 
asynchronism. 
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For healthy patients, the impedance vaJue at the beginning of QRS complex 
at ECG (Q-point), Z, is almost the same as the impedance value at the S-potnt, 2,. 
This impedance, difference is measured as Z^, (see FIG. 13). However, for ill 
patients, Z,., can be significant, (see FIG.14). Where Z„ is small, Kubicek's 
5 estimation of AZ=(dZ/dt) 1MX *VBT, is fairly accurate. For seriously ill patients 
exhibiting left-right ventricular asynchronism, however, the Kubicek equation 
generally underestimates AZ. Thus, the premise of the "prewave" should be 
accounted for in making the calculation. This invention compensates for the 
inherent underestimation of AZ in ill patients by adding the bioimpedance difference 
between S and Q points, Z t< to Kubkek's AZ estimation. By making this 
compensation, AZ is more accurately estimated, and the regression between 
bioimpedance and thermodilution cardiac output values is linearized. Thus, the 
computer system estimates AZ as 



10 



AZ = (dZ I dfl • ELVET ♦ Z, . 



15 where is the bioimpedance difference between S and Q points. 

After the computer system arranges all check points, it eliminates 
cardiocycles with certain aberrations. Fuzzy logic and fitness algorithms may be 
employed in this procedure. Several criteria are used for this purpose. The 
computer system first confirms that the time distances between the points described 

20 above for each cardiocycle (e.g. , point A, point T, Point S) do not exceed certain 
bounds. It also verifies that the amplitude difference between the start and the end 
of the cardiocycle do not exceed a predefined value. The amplitude change between 
the start and end of the cardiocycle must not exceed a predefined percentage of the 
maximum amplitude over that cardiocycle Furthermore, the ratio of the time- 

25 derivative bioimpedance signal amplitude at point A to the time-derivative 

bioimpedance signal at point T must be greater than a predetermined value; that is: 



mim) > c, 
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where Y(A) and Y(T) are the values of the. time-derivative bioimpedance signal at 
points A and T respectively, and "c" is the predetermined value. All cardiocycles 
that pass this stage are considered as "not very bad." The computer system then 
checks for the "neighbors criterion" to eliminate the effect of random noise in the 
bioimpedance signal. A three-dimensional "nearness" vector with elements (Ali, 
A2i, A3i) is calculated for each pair of cardiocycles in a 10 second time block. The 
individual elements of the nearness vector are determined using the following 
equations: 
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and 



where Y(x) is time-differentiated bioimpedance at a given point x, A, T, and B are 
check point positions in each cardiocycle (see FIG. 9); ST is the time between 
15 points S and T; and i and j are different cardiocycles. The computer system 
compares all of the calculated nearness vectors and eliminates those cardiocycles 
wherein the amplitudes of the nearness vectors exceeds certain thresholds. The 
comparison is made using a two-threshold analysts. If the nearness vector 
amplitude of two cardiocycles is less than a first predefined value, LI, the similarity 
is considered "good." If the similarity fails to be considered "good" (the similarity 
exceeds LJ), but is still less than second predefined value, L2, it is considered as 
"acceptable.* If the similarity fails to be considered "acceptable" for the 10 second 
block of data, the computer system compares the tested cardiocycle with up to 50 
previous "not very bad" cycles. If there is still no cycle similar to the one being 
considered, the cycle under consideration is regarded as noisy, and is rejected. If 
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the amount of -good" cardiocycles is big enough, all "acceptable" cardiocycles also 
rejected from (he final calculation. This methodology increases the stability of the 
calculations at the highest noise levels, as only "good" cardiocycles are used in the 
final calculations. 

5 After the bioimpedance signal processing is complete, the computer system 

performs the final calculation of hemodynamic parameters, together with means and 
variances for the entire data block. Heart Rate (HR) and Heart Stroke Volume (SV) 
are recalculated using the respective methods described above, except that the 
processed bioimpedance signal is used instead of the unprocessed bioimpedance 
signal. Cardiac Output (CO) is calculated as the product of HR and SV; that is: 

CO-SV-B1L 
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An outline of the general methodology of the pmsem ^ set foltfa ^ 

FIG. 6A. FIG. 6B depicts the flow chart of signal processing in the present 
invention. 

While the present invention has been described in terms of a preferred 
embodiment, those of ordinary skill in the art will recognize and appreciate that it is 
not so limited. Many additions, deletions and modifications to the disclosed 
embodiment may be made without departing from the scope of the invention as 
hereinafter claimed. 
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CLAIMS 

What is claimed is ; 

1 An apparatus for determining a subject's heart rate, heart stroke 
volume, and cardiac output from detected thoracic btoimpedance signals and 
5 electrocardiogram comprising an electrode array for detecting said thoracic 
bioimpedance signals and said electrocardiogram. 

2. The apparatus of claim 1 further comprising: 
an upper influencing electrode placed on the subject's head; 

10 a lower influencing electrode placed on the left lower extremity of the subject; 
an upper pair of detecting electrodes placed on the subject's neck; and 
a lower pair of detecting electrodes placed on the trunk of the subject. 

3. The apparatus of claim 2, wherein the placement geometry further 
IS comprises: 

an upper influencing electrode placed on the subject's forehead; 
a lower influencing electrode placed in the general area of the subject's left knee; 
a pair of upper detecting electrodes placed on the subject's neck; and 
a pair of lower detecting electrodes placed laterally on opposite sides of the 
20 subject's chest. 

4. The apparatus of claim 3, wherein said upper influencing electrode 
comprises a spot electrode for orientation on vertical and horizontal center lines of 
said subject's forehead. 

25 

5. The apparatus of claim 3, wherein said lower influencing electrode 
comprises a spot electrode, the placement of which satisfies the relationship L < 
5R, where L Is the vertical distance between said upper and said lower influencing 
electrodes and R is the radius of said subject's chest. 

30 

6. The apparatus of claim 3, wherein said upper detecting electrodes 
comprises a pair of spot electrodes oriented symmetrically on opposite sides of said 
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subjccfs neck .Hong a horizontal line approximately 4 centimeters above the base of 
said subject's neck. 



7 - The apparatus of claim 3, wherein said lower detecting electrodes 
5 further compri.se a pair of electrode assemblies, each assembly providing contact 
surface area between about 12 square centimeter* and about 30 square centimeters 
onented lateral* on opposite sides of said subject's chest at approximately xiphoid 
process level. 

10 8 - The apparatus of claim 7, wherein each assembly further comprises 

four spot electrodes, wherein each spot electrode comprises approximately 4 square 
centimeters contact surface area with each said spot electrode centered at comers of 
a square with sides measuring 5 centimeters, and all 4 spot electrodes of said 
assembly electrically connected to each other. 

15 

9. The apparatus of claim 8, wherein the top spot electrodes of each 
assembly lie on the xiphoid process level of the subject. 

10. A method for processing a bioimpedance signal and 
electrocardiogram for deriving heart rate, heart stroke volume, and cardiac output 
comprising: 

registering gain-phase-frequency (GPF) characteristics of input analog devices for 

measuring bioimpedance; 
registering gain-phase-frequency (GPF) characteristics of input analog devices for 
25 measuring electrocardiogram; 

measuring bioimpedance as a function of time over a given time period with said 
bioimpedance input analog devices and generating a bioimpedance goal; 
measuring electrocardiogram as a function of time over said given time period with 
said electrocardiogram input analog devices and generating an 
30 electrocardiogram (ECG) signal; 

correcting the bioimpedance signal for distortions based on the GPF characteristics 
previously registered; 
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correcting the electrocardiogram signal for distortions based on the GPF 
characteristics previously registered; 

determining valid QRS complexes associated with each cardiocycle of said 
electrocardiogram signal in the given time period; 
5 locating check-points on said valid QRS complexes; 

processing one of said ECG signal and said corrected bioimpedance signal to 
estimate heart rate; 

time-differentiating the corrected bioimpedance signal; 

determining check-points for each said cardiocycle of the time-differentiated 
10 bioimpedance signal in the given time period; 

determining effective left ventricular ejection time (ELVET) using said time 

differentiated bioimpedance check-points in relation to said corresponding 
QRS check-points; 

determining a novel correction factor Z M using said time differentiated 
15 bioimpedance check-points in relation to said corresponding QRS check- 

points; 

calculating stroke volume as a function of said ELVET, maximum time- 
differentiated bioimpedance (dZ/dt)^, specific blood resistivity (P), distance 
(L) between two bioimpedance voltage sensing electrodes of the 
20 biotmpetiance analog input device, baseline bioimpedance (Zq), said 

correction factor Z M , and a novel scale factor (K); and 

calculating cardiac output by multiplying said stroke volume by said heart rate. 

1 1 . The method of claim 10, wherein said registering gain-phase- 
25 frequency (GPF) characteristics of said bioimpedance input analog devices 

comprises determining phase- frequency and gain-frequency characteristics of a 
transducer employed in detection of said bioimpedance prior to use thereof in said 
detection. 

30 12. The method of claim 10, wherein said correcting of measured 

bioimpedance signal comprises digitally filtering and phase correcting said measured 
bioimpedance to remove distortion in the output of said transducer. 



WO 97/37591 



PCT/SG97/00013 



-28- 

13. The method of claim 10, wherein said registering gain-phase- 
frequency (GPF) characteristics of said electrocardiogram input analog devices 
comprises dctcimining pfiasc-fircquency and gain-frequency characteristics of a 
transducer employed in detection of said electrocardiogram prior to use thereof ii 
5 said detection. 



10 



15 



14. The method of claim 10, wherein said correcting of measured 
electrocardiogram signal comprises digitally filtering and phase correcting sai* 
measured electrocardiogram to remove distortion in the output of jaid transducer. 

15. The method of claim 10, wherein said estimating heart rate comprises 
using a power spectrum of the bioimpedance signal, and an auto-convolution 
function of the said power spectrum. 

16. The method of claim 10, wherein said estimating heart rate comprises 
processing the electrocardiogram signal. 



17. The method of claim 12, wherein said correcting of measured 
bioimpedance signal further consists of suppressing breath waves co remove 
20 undesired power spectra components and generate a bioimpedance signal of restored 
shape. 



18. The method of claim 10, wherein said deteimining valid QRS 
complexes comprises determining a distribution of all peaks measured in an 
electrocardiogram over a period of time (from the front (E,) and back (EJ 
amplitudes of the peaks, calculating the (E,, amplitude envelope and rejecting all 
peaks outside the envelope. 

19. The method of claim 11, wherein said determining transducer phase- 
frequency and gain-frequency characteristics comprises; 

generating a high precision sinusoidal impedance signal with peak-to-pcak 
impedance of approximately 0.2 Ohms and baseline impedance of 
approximately 100 Ohms to 200 Ohms; 
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connecting said sinusoidal impedance signal to said transducer; 
measuring the output from said transducer; and 

calculating a g;iin-phase-frequency characteristic, H(f), of the transducer in a 
predefined frequency range. 

5 

20. The method of claim 19, further comprising generating said high 
precision sinusoidal impedance signal using a voltage-to-irapedance converter 
including a photoresistor, a photoemitter, a power supply and an analog-to-digital-to 
analog computer interface. 

10 

21. The method of claim 20, further comprising producing through said 
interface a set of test signals with predefined amplitude and frequency range, 
detecting an output signal from said transducer, and analyzing said output signal to 
determine said phase-frequency and said gain-frequency characteristics. 

15 

22. The method of claim 19, further comprising employing posterior 
signal processing to correct linear gain-pbase-frequency distortions by converting 
real operating characteristics of the transducer to predefined characteristics, wherein 
phase shift is zeroed and gain is assumed to be constant in a predefined frequency 

20 range. 

23. A method of heart rate estimation, comprising: 
calculation of a power spectrum of a bioimpedance signal; 

multiplication of said power spectrum by a selected amplitude-frequency function to 
25 differentiate the signal and suppress breath harmonics; 

autoconvoluting the resulting power spectrum according to the formula 

ASW-PSoV)PSaQfrPSa&fl ... ; 



30 



and 

determining a maximum amplitude value of autoconvolution in a predefined 
frequency range as an estimation of heart rate. 
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24. A method for determining cardiocycles, comprising: 
filtering a bioimpedance signal to emphasize fronts of cardiocycles; 
calculating a time-amplitude envelope of said caidiocycles by analyzing the first five 

harmonics of the power spectrum of said bioimpedance signal after said 
5 filtration; 

selecting said cardiocycle fronts by comparison with said calculated time-amplitude 

envelope; and 
rejecting erroneously-detected fronts. 

10 25. A method of selecting valid cardiocycles from corrected 

bioimpedance signals to eliminate cardiocycles having interference artifacts, 
comprising: 

detecting time and amplitude relations referencing check points within individuals of 

a plurality of cardiocycles; 
IS comparing said time and amplitude relations between individuals of a said plurality 

of cardiocycles; and 
further examining selected cardiocycles which exhibit the presence of artifacts 

according to a plurality of comparison criteria. 

20 26. The method of claim 25, further comprising: 

constructing a multi dimensional vector for each selected cardiocycle; 
comparing said inulti-dimensional vector with such vectors for other cardiocycles 
and; 

rejecting the cardiocycles with vectors having no neighboring vectors when 
25 companxl to last 50 valid cardiocycles and other candidate cardiocycles. 



27, A method of deriving effective left ventricular ejection time from 
measured bioimpedance signal and measured electrocardiogram signal, comprising: 
filtering said measured bioimpedance signal and suppressing breath waves therein; 
30 filtering said measured electrocardiogram signal; 
detecting a valid cardiocycle; 

calculating the time-derivative of said bioinipedance signal Y(x); 
determining the maximum value of the time-derivative (dZ/dt)^; 
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determining effective ejection start time (S-point); 
determining effective ejection end time (T-point); and 
calculating effective left ventricular ejection time (EL VET) as change in time 
between effective ejection start time and end time. 

5 

28. The method of claim 27, wherein determining effective ejection start 
time comprises: 

determining the global maximum for a given valid caidiocycle of a time- 
differentiated bioimpedance signal and designating said maximum as point A; 
10 tracing back in time from corresponding point A on electrocardiogram to point S.; 
looking for abnormalities in the bioimpedance signal between points A and S a ; 
if there are no said abnormalities, then the cardiocycle is rejected as noisy; 
if there are any said abnormalities the one closest to point S b approached from the 
right is selected as the ejection start time S; 
15 otherwise the abnormality nearest S b approached from the left is selected as the 
ejection start time S. 



29. The method of claim 28, wherein determining effective ejection end 
time comprises: 

determining the first (T|) and second (T ? ) local minimums at the time-differentiated 

bioimpedance signal after point A; 
analyzing the depth of the signal curve at each of the first (T,) and second fTJ local 

minimums; 

if the depth of the second (T,) minimum is greater than a predetermined fractional 
value of the depth of the first (T,) minimum, selecting the second minimum 
(T 2 ) as T 0 ; 

otherwise, selecting T, as T 0 ; and 
identifying the T-point as the nearest local minimum before point T 0 on the graph of 
the curve generated by the second derivative of Y(x). 

30. The method of claim 29, further comprising, after identifying T, and 
T 3 and before identifying the T-point and regardless of relative amplitudes of T, and 
T,: 
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detecting a back edge of a T-wave in the electron cardiogram signal; 
determining if T, or T, is out of bounds of the back of the T-wave; and 
if so, selecting the one of T, or T 7 which is not out of bounds as T 0 . 

5 31. In the method of claim 28, abnormalities arc selected from the group 

comprising: dZ/dt zero crossing (Q-point), local rrwiirmim in dZ/dt, and local 
maximum in the third time-derivative of the bioimpedance signal, d 3 2/de. 

32. A method of determining stroke volume for a patient, comprising: 
10 determining specific blood resistivity P; 

measuring a distance L between two bioimpedance electrodes apphed to the patient; 
determining the base thoracic impedance Zoi 
determining EL VET; 

determining AZ, impedance changes due to blood influx; 
and calculating stroke volume SV according to the equation 



IS 



20 



25 



where K is a novel scale factor related to body composition of the patient. 

33. The method of claim 32, further comprising calculating K as 

K = ^ - K t • (SCBESTf(H*t • W^)\ 

where 

SCHEST - (POfESTL+PNECKPCHEST+PNEClPm*. 

34. Hie method of claim 33, wherein K,,, K„ K 2) K, arc gender and age 
dependent and lie in ranges of 
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JCo € [1-4]; K x € [3-16]; c [0-11; ^ e [0.1-2], 

35. A. method of claim 32, wherein determining AZ comprises; 
locating the start of QRS complex of an ECG signal and labeling it point Q; 
determining the impedance at point S, Z 4 ; 
determining the impedance at point Q, Z q ; 
5 calculating the impedance difference Z,. q between points S and Q; 
estimating AZ according to the formula 



LZ^idZim^ <ELVET + Z g _ q 



36, A method of breath wave suppression for a bioimpedance signal, 
comprising; 

10 calculating the Fourier transform of the signal; 

locating the first and second frequency harmonics of cardiocycles in the calculated 

spectrum of the signal; 
estimating the width of each of the harmonics; 

suppressing frequency harmonics below the lower bound of the second harmonic 
IS except for harmonics within the bounds of the first frequency harmonic; and 

calculating the inverted Fourier transform of the signal. 
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